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ABSTRACT 

We study the evolution of black holes (BHs) on the Mbh — cr and Mbh — Mb u i ge planes as a function 
of time in disk galaxies undergoing mergers. We begin the simulations with the progenitor black hole 
masses being initially below (A log A/eh, i ~ —2), on (A log A/eh, i ~ 0) and above (A log Mbh,! ~ 0.5) 
the observed local relations. The final relations are rapidly established after the final coalescence of 
the galaxies and their BHs. Progenitors with low initial gas fractions (/ gas = 0.2) starting below 
the relations evolve onto the relations (AlogA/BH.f ~ —0.18), progenitors on the relations stay there 
(A log Mbh, f ~ 0) and finally progenitors above the relations evolve towards the relations, but still 
remaining above them (A log Mbh, f ~ 0.35). Mergers in which the progenitors have high initial gas 
fractions (/ gas = 0.8) evolve above the relations in all cases (A log Mbh, f ~ 0.5). We find that the 
initial gas fraction is the prime source of scatter in the observed relations, dominating over the scatter 
arising from the evolutionary stage of the merger remnants. The fact that BHs starting above the 
relations do not evolve onto the relations, indicates that our simulations rule out the scenario in which 
overmassive BHs evolve onto the relations through gas-rich mergers. By implication our simulations 
thus disfavor the picture in which supermassive BHs develop significantly before their parent bulges. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: interactions — methods: numer- 
ical 



1. INTRODUCTION 

Observations in recent years have revealed a strong cor- 
relation in the local Universe between the central super- 
massive black holes (BHs) and their host galaxies as man- 
ifested in the relation between the B H mass and the bulge 
veloc i ty dispersion, Mbh ~ cr (e.g. [Ferrarese fc Merrittl 
120001 iTremaine et all l2002t iGiiltekin et all 120091: L the 
bulge stellar mass Mbh — Mh„io- (e.g. iMagorrian et "all 



the conc entration of light in 



1998; 



the galaxy (e.g. iGraham et al.|[200ll ) and the bulge bind- 
ing e nergy, Mbh — MbuigeC 2 (e.g. lAller fe Richstonel 
2007). The evolution of these relations with redshift is 
still unclea r with some stud i es fin ding evolution in the 
M RH - a (IWoo et all 120061 I 2008D and M RH - -M bllle r 
(|Peng et all l2006t iTreu et all 120071 : iDecarli et all 2009) 



relations, with the high redshift BHs being overmas- 
sive for a fixed a and Mb u i g c compared to the local 
relations, whereas other studies are consistent with no 
redsh i ft evolution in the observed correlations (|Gaskelll 
l2009t iJahnke et al.l [20091 ). A possible explanation for 
this discrepancy lies in the uncertainties in observational 
selection biases and in the evolution in the intrinsic 
scatter that is typ ically stronger for larger BH masses 
(|Lauer et al.ll2007h . 

The observed correlations are typically explained 
using theoretical models relying on some form of 
self-regulated BH mass growth, in which gas is fed 
to the central black hole until the black hole re- 
leases sufficient energy to unbind the gas and blow 
it away in mome n tum- or pressure-driven wind s (e.g . 
Silk fc Reed HoM iFabianl 11993 iBurkert fc Silld 120011 



Ciotti fc Ostriker^OOTHCiotti et al.ll2009D . The observed 



correlations and their evolution with redshif t have been 
reproduced in semi-analytic models (e.g. ICroton et al.1 
2006|; ISomerville et al.l |2008[ L in self-consistent numer- 



ical simulations of both isolate d galaxies and galaxy 
mergers (iDi Matteo et~all 120051: ISpringel etall l2005bi : 
iRobertson et al.l l2006f ) as well as in galaxies simu- 
lated in a full cosmological setting (Siiacki et al.ll2007l 
IDi Matteo etall 120081 : iBooth fc Schavdl2009h . The key 
assumption in these models is that the galaxies un- 
dergo a brief radiatively-efficient q uasar phase triggered 
by gas-rich galaxy merge rs (e.g. iTacconi et al.1 120021 : 
iHopkins et alll2006l l2008f) during which the bulk of the 
BH growth is taking place and the observed correlations 

are established. 

In a previous paper ([Johansson et al.l 12009. hereafter 
J09), we showed that the merger remnants of both equal- 
and unequal-mass mergers of disk and elliptical galaxies 
satisfy the observed Mbh — cr and Mbh — Mbuige cor- 
relations. In this Letter, we study for the first time in 
detail the evolution of the BH scaling relations during 
disk galaxy mergers using a new sample of high resolu- 
tion simulations including a self-consistent formulation 
for BH feedback. We seed the BHs initially with masses 
corresponding to locations below, on and above the ob- 
served relations. Thus, in contrast to previous studies 
we study here for the first time also overmassive BHs ly- 
ing initially above the observed relations. The BH scal- 
ing relations are defined for merger remnants that have 
reached their final dynamical state and it is not obvious if 
the scaling relations are valid during the merging process 
and at what stage the galaxies evolve onto the relations. 

2. SIMULATIONS 

The simulations were performed using the TreeSPH- 
code GADGET-2 (|Springell [20051) on the local Altix 
3700 Bx2 machine. The code includes radiative cool- 
ing for a primordial composition of hydrogen and he- 
lium. Star formation and the associated supernova 
feedback is included using the sub-resolution model of 
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TABLE 1 

The simulated merger sample 



Model 


MBH.il" 


MBH.i2 n 


Mass ratio 


a 


/gas 




-/Vdisk,tot 


N bl 


l.tot 


^DM.tot 


M BH / 


Tbul.f ° 


M bul , f d 


11B2BH 


1.0 


1.0 


1:1 


25 


0.2 


120 000 


480 000 


200 


000 


800 


000 


459 


183.8 


5.07 


31B2BH 


1.0 


1.0 


3:1 


25 


0.2 


80 000 


320 000 


133 


333 


533 


333 


271 


151.8 


1.90 


31B8BH 


1.0 


1.0 


3:1 


25 


0.8 


320 000 


80 000 


133 


333 


533 


333 


1136 


181.6 


2.97 


31B2BH1 


1.0 


1.0 


3:1 


100 


0.2 


80 000 


320 000 


133 


333 


533 


333 


82.5 


144.7 


1.26 


1102BH 


159 


159 


1:1 


25 


0.2 


120 000 


480 000 


200 


000 


800 


000 


600 


174.6 


3.38 


3102BH 


159 


36.4 


3:1 


25 


0.2 


80 000 


320 000 


133 


333 


533 


333 


253 


132.0 


1.55 


3108BH 


159 


36.4 


3:1 


25 


0.8 


320 000 


80 000 


133 


333 


533 


333 


1698 


155.2 


3.02 


11A2BH 


477 


477 


1:1 


25 


0.2 


120 000 


480 000 


200 


000 


800 


000 


1200 


164.2 


3.35 


31A2BH 


477 


109 


3:1 


25 


0.2 


80 000 


320 000 


133 


333 


533 


333 


652 


128.7 


1.59 


31A8BH 


477 


109 


3:1 


25 


0.8 


320 000 


80 000 


133 


333 


533 


333 


842 


122.3 


1.21 



"Initial BH mass in 1O 5 M 

6 Final BH mass in 10 5 M Q 

c Final stellar velocity dispersion in km/s 

d Final bulge mass in 1O 1O M0 




Fig. 1. — The relative distance between the BHs (top), the total 
SFR (middle), and the evolution of the total BH accretion rate 
(bottom) as a function of time. The star symbol indicate the time 
of BH coalescence and the filled circles at the bottom of the plots 
indi cate the time at which (AfBH^-Mbulge) are evaluated in Figs. 

EH 



ISpringel fc Hernquistl ([2003). I n this model the ISM is 
treated as a two-phase medium (jMcKee fc Ostrikerlll977t 
I Johansson fc Efstathioul |2006[ ) in which cold clouds are 
embedded in a tenuous hot gas a t pressure equi l ibrium . 
BH feedback is modeled using the ISpringel et alj (|2005bl ) 
effective model in which the BH sink particles accrete 
gas from the surrounding medium according to a Bondi- 
Hoyle-Lyttleton parametrization 



Mb = 



( c 2 +u 2)3/2 ' 



(1) 



with the maximum accretion rate set by the Eddington 
limit. Here p and c s are the density and sound speed 
of the surrounding gas, respectively, v is the velocity of 
the BH relative to the surrounding gas and a is a di- 
mensionless efficiency parameter. Following J09 a total 
of 0.5% of accreted the rest mass energy is then injected 
as thermal energy into the gas surrounding the BH parti- 
cle. Throughout this paper we use the BH repositioning 
method that ensures rapid merging of the BHs once the 
two BHs reside in the same parent galaxy. We assume 
that the BHs merge instantly if they enter the smoothing 
length of each other and if their relative velocity is below 
the local sound speed. 

The progenitor g alaxies are setup following the method 
of ISpringel et al.l (|2005bh and as detailed in J09 with 
all primary galaxies having v v - 1T — 160 kms^Vvir = 
160 /i _1 kpc corresponding to a virial mass of M v - lr — 
9.53 x lO 11 /^ 1 ^?© 1 analogous to the Milky Way. The 
dark matter profiles i n all model galaxies are described 
by a iHernquistl (|1990f) profile constructed using the con- 
centr ation parameter c = 9 of the corresponding NFW 
halo H avarro et al.lll997t ). The dark matter halos are 
populated with exponential disks with a baryonic mass 
fraction of irid — 0.041 and a stellar Hernquist bulge with 
a fraction mi, — 0.01367 (mj = 1/3™^) of the total virial 
mass, M v ; r . The disk is setup with a fractional gas con- 
tent of /gas, with the rest being disk stars, with the scale 
lengths of the disk and bulge computed as detailed in 
J09. The secondary models in the unequal-mass mergers 
are identical to the primary models, except for the fact 
that all components are scaled down in mass by a factor 
of three. 

In setting up the primary galaxy disks we use a rela- 
tively high numerical resolution of 300,000 disk particles 
of which a fraction / gas are gaseous, 100,000 bulge parti- 
cles and 400,000 dark matter particles. Each gas particle 
can only spawn a single stellar particle resulting in a con- 
stant particle mass resolution of m Da r = 1.30x 10 5 h~ 1 M & 
for all baryonic particles and ttt-dm = 2.25 x 10 6 H~ 1 Mq 
for the dark matter particles. The gravitational softening 
length was set to e = 0.02/i _1 kpc for all baryonic parti- 
cles and the BHs and to e = 0.083/i _1 kpc for the dark 
matter particles. Finally we insert a seed black hole at 
rest in the center of each galaxy model, where we vary 
the BH mass depending on if the models should lie ini- 

1 h = 0.71 is defined such that H = lOO/i kms _1 Mpc _1 
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Fig. 2. — The evolution of the BHs in the B-merger sample ( A log Mbh. i ~ —2) as a function of time on the Mbh — cr plane (top left 
panels) overplotted by lines giving the observed correlation from ITremaine et a l, (2002). The bottom left panel gives the evolution on the 
Mbh — Mb u ige plane overplotted by the observed correlation from lHaring fc R ix (2004). Each point is separated by At = 0.2 Gyr, with 
filled circles indicating the primary galaxy, triangles indicating the secondary galaxy, stars showing the time of BH coalescence and filled 
squares the final state. The panels on the right give the logarithmic offset (Eq. [2} from the corresponding observed relation as a function 
of time. 



tially below (B-series), on (O-series) or above (A-series) 

the observed local BH scaling relations. 

We adopt orbital geometry G13 (|Naab fc Burkertl 
2003) for all simulations in this paper. This geome- 
try corresponds to the inclinations i p = — 109, i s = 180 
and the arguments of the pericenter tu p — 60, uj s = 
for the primary and secondary galaxies, respectively. 
The galaxies approach each other on a parabolic or- 
bit where the initial separation of the progenitors is 
-Rinit = 0.5(r v i r! p + r v i r!S ) and the pericentric distance is 
J"peri = r d ^ P + r diS , where r viriP , r d!P and r vir , s , r d!S are the 
virial and disk scale radii for the primary and secondary 
galaxies, respectively. All simulations were evolved for a 
total of t = 3 Gyr with the merger typically taking place 
at t ~ 1.5 Gyr. The main numerical parameters and final 
properties of all 10 simulations are summarized in Table 

The parameters governing the multi-phase feedback 
are as detailed in J09 resulting in a star formation rate 
(SFR) of ~ lMoyr -1 for our primary galaxy. The ac- 
cretion parameter a (Eq. [T]) can essentially be seen as 
an empirical correction factor that translates from the 
resolvable low mean density to the time-averaged small- 
scale density at the location of the BH. We found that 
for a = 25 the BHs reach Eddington limited mass growth 
during the merger. This is not the case for the value 
a = 100 used in J09 (see Fig. [TJ. The correspond- 
ing evolution on the Mbh — c plane is such that the 
3B2BH1 (a = 100) simulation does not evolve onto the 
relation, whereas the 3B2BH (a = 25) evolves onto the 
relations reproducing the lower resolution results from 
J09 (see Fig. [2]). We conclude that as the spatial reso- 



lution in the simulations is increased the corresponding 
a value needs to be decreased in order to reproduce the 
lower resolution runs. Thus, we adopt a = 25 for all 
subsequent simulations. Finally, we note that our results 
do not depend critically on a as long as the condition 
of Eddington limited mass growth during the merger is 
fulfilled. 

3. METHODOLOGY 

We plot in Fig. [T] the relative distance between the 
BHs, the total SFR and the total BH accretion rate for 
a subsample of six simulations. In addition, we overplot 
filled circles at time intervals of At = 0.2 Gyr through- 
out the 3 Gyr simulation indicating the time at which 
the galaxies are extracted for analysis. This gives a to- 
tal of 16 snapshots during the merger. We perform the 
analysis using all stellar particles found in a sphere with 
a radius of r = 30 kpc centered on the BH of the corre- 
sponding galaxy. The repositioning technique employed 
in this study ensures that the BHs always trace the high- 
est density central structure of their parent halo. 

The inclusion of a massive BH initially has a marked 
effect on the star formation and BH accretion histories 
of the mergers (see Fig. [IJ. The higher the initial BH 
masses the lower the resulting SFRs are throughout the 
merger. The SFRs are suppressed by a factor of ~ 4 
in the O-series and a factor of ~ 8 in the A-series sim- 
ulations compared to the B-merger series. In addition, 
massive BHs suppress the peak BH accretion rate during 
final coalescence by a factor of > 10. 

Following observational estimates we calculate from 
each snapshot the line-of-sight stellar velocity dispersion 
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Fig. 3. — The evolution of the O-merger sample with the BHs starting initially on the local scaling relations (A log AfeH,i ~ 0). Symbols 
and the observed relations as denned in Fig. [2] 



a inside the effective radius using 50 randomly projected 
realizations of the galaxies. Together with the extracted 
BH masses we then plot the evolution of the BHs during 
the mergers on the Mbh — & plane (top left panels in 
Figs. EE]). Here we also overplot the observed local rela- 
tion (log(M BH /M ) = a + frlog (cr/cr n ), a = 8.13 ± 06, 
b = 4.02±0.32, cr = 200 kms- MTremaine et"aT]|2002D as 
a solid line with the dashed lines indicating the la errors. 
In the top right-hand panels of Figs. [50 we calculate as 
a function of time the relative offset of the simulated BH 
masses from the ones predicted by the relation defined 
as, 

A log M BH = log Mbh, simulation 

- log M B h, predicted (2) 

where a positive A log Mbh indicates an overmassive and 
a negative A log Mbh an undermassive BH with respect 
to the observed correlation (the Mbh — a relation in this 
instance). The dashed lines in the plots give the mean 
la error in the velocity dispersion range probed by our 
simulations. 

In the bottom panels of Figs. EE1 we plot the evolu- 
tion of the BHs on the Mbh — Mb u i gc plane, with the to- 
tal bulge masses derived following the method outlined 
in iNaab fc Trujillol (]2006h . Artificial images smoothed 
with a Gaussian filter were created of every snapshot 
seen from 50 random projections by binning the cen- 
tral 30 kpc into 128x128 pixels. We then fit both a 

pure Sersic (S(r) = X(0)e~ b "( r / rc ' )1/ ' 1 ) as well as an ex- 
ponential disk and Sersic bulge component simultane- 
ously (E(r) = £ I) (0)e r /' lD +I] s (0)e- b "(''/ r ^ 1/ " f3 ). Here 
r e ,r e s are the effective radii and n, the Sersic in- 
dices of the bulge component, Hd is the scale length of 
the disk, E(0) are the central surface densities and b n 
is chosen so that r e encloses half the total mass. We 
then solve for the best fit bulge-to-total (B/T) ratio for 



each projection and give the total stellar bulge mass as 
the mean of the best fit values. In addition, we over- 
plot the observed local Mbh — Mbui go relation with errors 
(log (Mbh /Mq) = c+dlog(MJ10 11 M w ), c = 8.20±0.10, 
d = 1.12 ± 0.06, lHaring fc R ix 2004) together with the 
corresponding evolution of A log A/bh as a function of 
time in the bottom right panels. 

4. RESULTS 

4.1. Mergers with A/bh starting below the relations 

We plot in Fig. Elthe evolution of our B-merger sample 
on the planes of the two BH scaling relations. For this 
sample the initial BH mass is always set to Mbh,! = 
10 5 Mq (A log M B h, i ~ -2) similarly to J09 and thus 
rendering the BH initially dynamically unimportant. 

All models show initially relatively modest evolution 
towards the observed relations, with the first passage at 
t ~ 0.6 Gyr primarily only seen as a mild increase in a 
and Mbuigo- In the phase t ~ 0.6 — 1.2 Gyr when the 
galaxies are falling back towards the second encounter 
both a and Mb u i go are increasing more rapidly, due to 
stronger merger- induced tidal disturbances. After the 
second encounter the galaxies rapidly merge within t ~ 
0.2 — 0.4 Gyr, with the shorter timescales corresponding 
to the equal-mass merger and the merger with the higher 
initial gas fraction of / gas = 0.8 (see Fig. [T]). 

This e volution is qualitat i vely similar to the ea rly re- 
sults of INaab et ail (|2006[) ; iDasvra et all (|2006[) . who 
found that the velocity dispersion increases primarily in 
two phases, with mild increase after the first encounter 
an d very rapid evoluti on at the time of the coalescence. 
In lDasvra et al.l (|2006f ) the authors also conjectured that 
the BHs would fall on the Mbh — o already after the first 
encounter if a constant fraction of 1% of the central gas 
supply would be accreted onto the BH. However, these 
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Fig. 4. — The evolution of the A-merger sample with the BHs starting initially above the local scaling relations (A log Mbh, i ~ 0.5). 
Symbols and the observed relations as defined in Fig. [2] 



early studies were purely hydrodynamical, lacking self- 
consistent BH growth and feedback. 

We find that although the BHs evolve towards the re- 
lation during the merger, the final relations are only es- 
tablished at the time of the final coalescence of the BHs. 
Typically the velocity dispersions of the remnants are 
somewhat too large immediately after the merger. The 
effective velocity dispersion then slowly decreases due to 
some low-level residual star formation with the stars pref- 
erentially settling in a disk-like structure, with this effect 
being stronger in high-gas fraction mergers. However, 
this effect is less strong in low-gas fraction and low mass 
ratio mergers, as star formation and BH mass growth 
are more strongly terminated in these mergers after fi- 
nal coalescence (Fig. [1] see also ISpring el et al.l l2005a: 
Uohansson et al.ll2008T) . On the other hand, we find that 
the bulge mass is established on a short timescale imme- 
diately after coalescence with the pre-coalescence bulge 
mass typically being markedly lower. 

The prime driver for the scatter is caused by the ini- 
tial gas fraction / gas of the progenitor galaxies, with the 
scatter produced by the evolutionary stage of the merger 
remnants after the final coalescence being only a sec- 
ondary effect, as can be seen in the relatively modest 
evolution between coalescence (star symbols) and the fi- 
nal merger remnant state (square symbols). 

4.2. Mergers with Mbh starting on the relations 

In Fig. [3] we study the evolution of the O-merger 
sample, for which we set the initial seed BH mass such 
that the model galaxies lie on the observed relations 

(A log Mbh, i ~ 0)- 

Both the primary and secondary galaxies in the merg- 
ers typically evolve initially along the observed relations, 
but after the first encounter before final coalescence the 



secondary galaxy typically evolves horizontally off the 
relations, again due to the tidal disturbances after the 
first passage (Fig. [3]). However, after the merging of 
the BHs all remnants again evolve rapidly onto the rela- 
tions, with the exception of the high gas fraction merger, 
which evolves off the relation. Thus for progenitors with 
initially low gas fractions (/ gas = 0.2) the BH feedback is 
able to self-regulate the growth of the BHs maintaining 
them on the relation. For very high gas fractions the BH 
growth is too efficient with the self-regulation being un- 
able to maintain the galaxies on the observed relations. 

4.3. Mergers with Mbh starting above the relations 

Finally, m otivated by the observations of I Woo et alj 
we run simulations in which the BHs were 
initially overmassive by a factor of three (A log Mbh ~ 
0.5). 

The evolution of this A-merger sample is shown in Fig. 
2J Again we see initially weak evolution with the ini- 
tial offset maintained, but after the first encounter the 
secondary BH evolves horizontally towards the observed 
relation and in some instances even below the relations. 
After the merger all remnants again evolve above the re- 
lations. However, the final merger remnants of the low 
gas fraction (/ gas = 0.2) simulations have a lower off- 
set (A log Mbh) by ~ 0.2 dex compared to the initial 
state. This effect is stronger in the 1:1 merger, which 
is experiencing a more violent encounter, with more gas 
channeled to the center and thus stronger self-regulation 
of the BH. The high gas fraction merger (/ gas = 0.8) on 
the other hand evolves even further away from the ob- 
served relations increasing the offset typically by ~ 0.2 
dex. 

5. DISCUSSION 
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In this Letter, we have studied the evolution of galaxy 
mergers on the Mbh — o and -Mbh — Mbuige planes as 
a function of time. We have shown that progenitors 
with low initial gas fractions (/ gas = 0.2) starting be- 
low the relations evolve onto the relations, progenitors 
on the relations stay there and progenitors above the re- 
lations evolve towards the relations, but still remaining 
above them. Progenitors with high initial gas fractions 
(/gas = 0.8) evolve above the relations in all cases, with 
the initial gas fraction thus being the p rime source of 
scatt er in the observed relations (see also lHopkins et all 
2007). The evolution for all mergers is initially slow with 
the observed relations typically being rapidly established 
during a relatively short phase centered at the time of fi- 
nal coalescence of the BHs. 

The observations of IWoo et all (|2006t 120081 ) indicate 
that the BHs at high redshifts were typically overmassive 
by a factor of ~ 3 for a fixed a and Afbuigc compared to 
the local relation. These BHs could plaus ibly have been 
forme d during a brief quasar phase (e.g. I Hopkins et all 
2006) triggered during the mergers of very gas-rich galax- 
ies, thus resulting in overmassive BH masses for their 
given a and -Mbuige, as seen in our / gas = 0.8 simulation 
series. However, this being the case it is not obvious how 
these galaxies would evolve onto the local observed rela- 
tions until the present day. Another binary merger with 
massive BHs in place does not bring the galaxies onto the 
relations (O- and A-series), with high gas fractions merg- 
ers moving them even further away from the relations. 
Potentially, the bul ge mass could be increased by dry 
minor mergers fe.g. iNaab et al.ll2009f ). whereas internal 
secular processes could be responsible for increasing the 
velocity dispersion. However, it is not obvious how both 
the velocity dispersion and the bulge mass could be in- 



creased at a fixed BH mass. Finally, another possibility 
could be that some of the BH s undergoing merg ers are 
ejected in a sling-shot effect ()Saslaw et all 11974) . thus 
decreasing the total BH mass. 

The fact that the BHs starting below the relations 
evolve onto the relation, whereas the ones above do not, 
indicates that our simulations rule out the scenario in 
which overmassive BHs evolve onto the relations through 
gas-rich mergers, whereas undermassive BHs can evolve 
onto the relations in galaxy mergers. Thus, given the nu- 
merical limitations inherit in our simulations, our results 
disfavor the picture in which supermassive BHs develop 
significantly before their parent bulges. 

Finally, our current BH accretion and feedback pre- 
scription seems to describe the growth of BH during the 
merger phase adequately. However, it is not obvious that 
the present description is also valid during the secularly 
driven low accretion phase after the merger is completed. 
Some initial steps to improve the prescription have been 
taken by developing a n entirely new mom entum driven 
feedback p rescription ([ DcBuhr ct al. 2009). Recent ob- 
servations dSchawinski et all 120071 12001 " have indicated 
that there is a time lag of ~ 0.5 Gyr between the peak 
of star formation and the onset of AGN activity. Our 
current model has difficulties in reproducing this and 
thus developing models that shed light on this discrep- 
ancy might also ultimately help us understanding in how, 
when and why the BHs in galaxies evolve onto the ob- 
served relations. 

The numerical simulations were performed on the lo- 
cal SGI-Altix 3700 Bx2, which was partly funded by the 
Cluster of Excellence: "Origin and Structure of the Uni- 
verse" . 
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